
From the depths of scientific mystery end a baste knowledge of 
mathematics, weird and wonderful designs can be produced, from a 
computer with graphical output. Barry Mart, n and Mtke Mudge present 
a kaleidoscope of patterns and figures. 



'Chaos' is a respectable scientific dis- 
cipline emerging in many areas, 
ranging from turbulent flow through 
population dynamics, heart cell dyna- 
mics and neurological patterns to the 
study of lasers. Any readers in- 
terested in the underlying theory can 
refer to the publication, The Univer- 
sality of Chaos , a collection of papers 
edited by Predrag Cvitanovic and 
published by Adam Hilger (1984). 
However, with no knowledge of the 
theory and only a minimal know- 
ledge of basic mathematics, patterns 
and designs of breathtaking beauty 
and intricacy can be produced — all 
that is required is a microcomputer 
with graphical output. Colour 
graphics, although not essential, can 
be used to produce subtle effects 



and attractive designs. 

The procedure is based on itera- 
tions of ordinary real numbers. A 
starting point (xo,yo) * s chosen and 
then successive points (x^yi), (x 2 ,y 2 h 
and so on, are generated using some 
mathematical function. A digital com- 
puter plots the points so quickly that 
it appears as if electronic raindrops 
are falling onto the graphics screen. 
Eventually, after a thousand or so 
points have been plotted, a pattern 
begins to emerge. Patterns can be 
varied by changing the starting point 
or by changing numerical parameters 
in the mathematical function, thus 
giving infinite variety from any one 
point generator. 

Do not, however, get the impress- 
ion that the formulae have to be 



complicated in order to produce exn 
tic patterns. Very simple formulae ut 
ing, for example, functions such 
SGN(x) and SQR(ABS)x)), which a 
available on all micros, can yield ex 
quisite designs as can be seen from 
the detailed examples which ann^r 
in this article. ppear 



Some patterns seem to build up j n 
one area of the screen for thousands 
of iterations and then, suddenly, the 
whole process jumps to another re- 
gion and a further segment of the 
pattern develops. Obviously, some 
experimentation with scaling j s 
usually necessary to reveal the full 
content of a given point generator 
Increasing the magnification of an 
area surrounding a particular point of 
a pattern will reveal more of the 







microstructure of the pattern and will 
probably produce a thing of beauty 
in itself. 

Colours may be used to enhance 
the artistic appeal of the patterns, for 
example, by assigning colour accord- 
ing to how far the point being plot- 
ted is from the origin of coordinates. 
Alternatively, you can change colour 
according to the number of iterations 
required to compute the position of 
the point — for example, vary the 
colour every 100 iterations. 

The techniques for developing the 
patterns described in this article were 
inspired by the Mandelbrot sets pre- 
sented by Dr MM Novak and Jack 
Weber in the article 'Fractal Sets', 
PCW, December 1986. Mandelbrot's 
patterns emerge from complex num- 
bers (not easy things for many peo- 
ple to manipulate on their compu- 
ters) and the seeds for the iterative 
process are the points, infinite in 
number, found in a region of the 
plane. For every point on the screen 
a complicated calculation needs to 
be done many times over, and the 
computer generation of a complete 
pattern may take hours or even days, 
although Weber does describe ways 
to reduce this computation time. 

Furthermore, to obtain detailed re- 
solution of the final picture, high- 
powered graphics facilities are re- 
quired. The process described here, 
on the other hand, uses only real 
numbers and the pattern grows from 
only a single seed, the starting point 
(x 0 , Vo)- The pattern begins to appear 
after a minute or so of computation, 
and exciting pictures can be pro- 
duced using nothing more elaborate 
than the screen of a portable TV set. 

Using such pattern-generating 
techniques, it is possible to design 
your own personalised wallpaper or 



textiles so that no-one else in the 
world would be likely to have the 
same pattern. Imagine going into a 
shop, keying a couple of your own 
secret numbers into a computer ter- 
minal, and a computer-based auto- 
matic 'pattern designer' then pro- 
duces rolls of unique personalised 
wallpaper or even a unique dress 
fabric. 

However, to return to the 'chaos' 
of reality upon which all of this is 
based, before you all branch out into 
computerised wallpaper design, we 
will look at a one-dimensional case 
(no pretty pictures) which will illus- 



trate how chaotic phenomena rise 
from an ordered situation. 

One-dimensional chaos 

As an example we will look at the 
equation (known to mathematicians 
as the logistic equation) (1) 
. . .x r +1=4Cx r (1-x r ), where r = 0,1,2 
. . . successively. We take a starting 
value x 0 anywhere between 0 and 1 
and then compute x 1# x 2 , x 3 . . . using 
equation (1). 

A procedure such as this is called 
an 'iterative scheme' and the com- 
putation of each successive x — 
value is called an iteration. Before 



10 


REM PROGRAM "hopalong". 




20 


REM THIS GENERATES PATTERNS BY ITERATION FROM "HOPALONG". 


40 






50 


MODEl 




60 


REM READ IN DATA: CONSTANTS P,Q,R, INITIAL POINT, CENTRE POSITION, 


70 


REM SCALE AND NUMBER OF ITERATIONS. 




80 


INPUT "CONSTANT P ",P 




90 


INPUT "CONSTANT Q ",Q 




100 


INPUT "CONSTANT R " ,R 




110 


INPUT "X COORDINATE OF STARTING POINT 


",X 


120 


INPUT "Y COORDINATE OF STARTING POINT 


" , Y 


130 


INPUT "X COORDINATES OF CENTRE POINT " 


,XC 


140 


INPUT "Y COORDINATES OF CENTRE POINT " 


,YC 


150 


INPUT "HORIZONTAL SCALE ",XS 




160 


INPUT "VERTICAL SCALE ",YS 




170 


INPUT "NUMBER OF ITERATIONS REQUIRED " 


>N 


180 


CLS 




190 


IC0UNT=0 




200 


T=1/SQR(2) 




210 






220 


REM START ITERATION. 




230 


FOR 1-1 TO N 




240 


X1®Y-SGN(X)*SQR(ABS (Q*X-R) ) 




250 


Y=P-X 




260 


X=X1 




2 70 


IC0UNT=IC0UNT+1 




280 


PRINT TAB ( 0,0) ; ICOUNT 




290 


REM ROTATE PLOT THROUGH -45 DEGREES. 




300 


U=T*X+T*Y 




310 


V=-T*X+T*Y 




320 


PLOT 69 ,XC+XS*U, YC-*-YS*V 




330 


NEXT I 




340 


STOP 
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graphics 



^m‘rEAD in data: CONSTANT A INITIAL POINT, 'CENTRE' POSITION, 

S “ale and number of iterations. 



1° ™ ITERATES WHICH "TILE" THE PLANE 

ll REM WITH AXEHEAD LIKE FORMS. 

50 

60 MODE1 
70 

80 KX-n jwnu- 

90 “SE "^COORDINATE ' OF STARTING POINT ",X 

100 INPU I „ COORDINATE of starting point ",y 
“S "X coordinates of centre point »,xc 

no llm "Y COORDINATES OF CENTRE POINT ,YC 
140 INPUT "HORIZONTAL SCALE ,XS 
I crv tmpmt "VERTICAL SCALE 
leo INPUT "NUMBER OF ITERATIONS REQUIRED 

170 CLS 
180 IC0UNT=0 
190 T=1/SQR(2) 

200 

210 REM START ITERATION. 

220 FOR 1=1 TO N 
230 Xl=Y-SIN(X) 

240 Y=A-X 
250 X=XI 

260 IC0UNT=IC0UNT+1 
270 PRINT TAB(0,0) ; ICOUNT 
280 REM ROTATE PLOT THROUGH -45 DEGREES. 

290 U=T*X+T*Y 
300 V=-T*X+T*Y 
310 PLOT 69 ,XC+XS*U,YC+YS*V 
320 NEXT I 
*nn stop 



starting, however, we must decide 
on a value for the constant C, as the 
behaviour of the iterative scheme is 
dependent on the constant we 
choose. Firstly, C must lie between 
and 1; this keeps the successive iter- 
ates within the range zero to one. 

C is less than 0.25, the values of x r as 
r becomes large approach zero. If C 
lies between 0.25 and 0.75, the 
values of x r approach 1-V4C. 

For C between 0.75 and 0.9, more 
interesting features begin to emerge. 
For example, when C=0.76, x r alter- 
nately approaches the two values of 
0.589356083 ... and 0.730591286 . . . 
and, as larger values of C less than 
0.9 (or, more accurately, 0.8925 . . .) 
are chosen, x r simultaneously 
approaches an increasing number of 
these 'limit' values, as shown in the 
box below. 



c 

0.2 
0.24 
0.76 
0.87 
0.89 
0.892 

Limit values 



The algebraic significance of this 
phenomenon is well-understood by 
the mathematical pundits and is, in- 
deed, quite predictable and orderly. 
However, as C gets closer to the spe- 
cial value 0.8925 . . ., order ceases 
and chaos takes over. The values of 
x r become distributed between 0 and 
1 in a thoroughly haphazard way, but 
order can return as C is further in- 
creased towards unity. For example, 
if you try X 0 =0.96, you should find 
that just three points are visited. 



All the various features described 
above can be displayed using the 
program 'CHAOSMW' on page 122, 
which is written in BBC Basic. This 
first draws the parabola having equa- 
tion (2) y=4Cx(1-x), together with 
the straight line having equation (3) 
y=x. It then carries out the iterative 
process which is illustrated in Fig 1 
where, in this case, C=0.8 and 
x o =0.3. 

CHAOSMW originated from 
Matthew Wells, who produced it as 
part of an undergraduate project at 
the University of Aston, Birmingham 



in 1986. However, he is making no 
claims to be the first to write such a 
program: many mathematicians dur- 
ing the last few years have con- 
structed similar routines. 

To help you centre and scale the 
plots on your own monitor, it will be 
useful to know that on a screen of 
1280 X 1024 pixels, xo = 200, yo = 
100, xs = 900 and ys = 900 gives a 
good picture. 

Two-dimensional chaos 

The two-dimensional mappings 
which we will consider are of the 



Number of 'limit' values 
1 
1 
2 

4 

8 

16 



form (4) x r +1=f(x r +y r ), y r +1=A-x r , 
where A is a numerical constant, r = 
0,1,2 . . . successively and starting 
values (x 0 , yo) are to be specified, 
and f(x,y) is a non-linear function of 
the ordered pair (x,y). Readers have 
the opportunity to investigate the 
consequences of replacing 
y r +1=A-x r by y r +1 =g(x r +y r ), a 
second non-linear function. Brief ex- 
perimentation within this area, 
however, suggests that the resulting 
patterns may be over-elaborate with 
a consequent loss of aesthetic 
appeal. 

The first family of mappings to be 
considered is (5) x r +1=y r -sin(x r ). 



GRAPHICS 



* o suC" 

The necessary program AXHEADS , 

X' « A-.-3 1415026 

.... this program displays repM 
tion of the four points (0,0), (0,77), 
Zrr) and M). This phenomenon is 
somewhat analogous to the predict- 
able and orderly behaviour previous- 
ly described for the one-dimensional 
mapping equation (1) when the con- 
stant C is less than 0.87 . . . We are 
not, therefore, surprised by the 
approximation to four-fold symmetry 
that results from choosing a value 
for A which is sufficiently chose to tt 



10 

20 

30 

40 

50 

60 



100 
110 INPUT 
120 INPUT 
130 INPUT 



REM PROGRAM "CHAOSMW". t^tcttp FDI1ATTON 

REM THIS PROGRAM DISPLAYS SUCCESIVE ITERATES OF THE LOGISTIC EQUATION 

REM X(N+1 )=4*C*X(N)*( l-X(N) ) . 

REM TO START IT ASKS FOR A VALUE OF C BETWEEN 0 AND 1 
REM AND FOR AN INITIAL VALUE OF X BETWEEN 0 AND 1. 

REM THIS PROGRAM WAS WRITTEN BY MATTHEW WELLS. 1986. 

70 MODE 1 
80 

90 INPUT "CONSTANT 0<C<1 ",C 

INPUT "STARTING VALUE 0<STARTX<1 ",STARTX 

MAXIMUM NUMBER OF ITERATIONS ",MAXIT 
X COORDINATE OF ORIGIN ",X0 
Y COORDINATE OF ORIGIN ",Y0 
140 INPUT "LENGTH IN PIXELS OF X-AXIS ",XS 
150 INPUT "LENGTH IN PIXELS OF Y-AXIS ",YS 
160 

170 REM DRAW AXES AND THE PARABOLA 4*C*X*(1-X). 

180 CLS 

190 PRINT TAB (0,0) 

200 PRINT TAB (4 , 31 ) "0" 

210 PRINT TAB(0 , 0) 

220 PRINT TAB (34 , 31 ) "1" 

230 PRINT TAB (0,0) 

240 PRINT "C=" ; C 
250 PRINT "X(0)=" ; STARTX 
260 PLOT 4 ,X0,Y0 
270 PLOT 5 , XO+XS , YO 
280 PLOT 4 ,XO,YO 
290 PLOT 5,X0,YS 
300 PLOT 4 ,X0, Y0 
310 PLOT 5 , XO+XS , YO+YS 
320 X=0.0 
330 SCALE=1/YS 
340 FOR 1=1 TO -YS 
350 X=I*SCALE 

360 PLOT 69 , X0+X*XS, YO+YS*4*X*C*( 1-X) 

370 NEXT I 
380 

390 M=0 
400 X=STARTX 

410 PLOT 4 , X0+STARTX*XS , Y0 
420 GCOL 0,1 

430 REM START ITERATION. 

440 Y=4*X*C*( 1-X) 

450 

460 REM PAUSE TO SEE DEVELOPMENT OF PLOT. 

470 FOR N=0 TO 400 
480 NEXT 

490 PLOT 5 , X0+X*XS , Y0+YS*Y 
500 PLOT 5 , XO+Y*XS , YO+Y*YS 
510 X=Y 
520 M*M+i 

530 REM INTRODUCE COLOUR CHANGE AFTER 50 ITERATIONS 
540 IF M<50 THEN GOTO 560 
550 GCOL 0,2 

560 IF M=MAXIT THEN STOP 
570 GOTO 440 



(see Fig 2 referring toX^Tu'hT 

20,000+ points have been D l 
and the colour changed eve 0tted 
iterations. However, the ^ ^ 
associated with the value Al\ Ults 
and using different numbers of > 15 
tions are displayed in Fig s 3 a '^ ra ' 
These pictures display some 0 f 4 ' 
infinite variety associated W j th ![j e 
point-generator equation (5) r as * he 
cussed in the introduction. * d,s ' 

Scaling with A=3.12, and in;*- , 
point. (0,0) on the screen mention^ 
above, is reasonable with y c ~R^ 
yc=550, xs=10 and ys=10 

A second family of mappings to be 
considered is (6) x r +1=y r - sinn ,;® 
x (Qx r -R|, y r+ 1=P- Xr , where 9 ?^ 
and R are numerical constants 
r=0,1,2 . . . successively and starting 
values (x 0 ,y 0 ) are to be specified 9 

The case when x 0 =0, y 0 =0, R~ n 
and Q=1 has been discussed by Mar 
tin. The Mathematical Gazette 
volume 70, number 452, June 1986 
pages 140-142, for a range of values 

The general investigation of the 
point generator, equation (6), may be 
carried out using the program 
'HOPALONG' on page 119 0 f this 
article (the name was coined by AK 
Dewdney who writes on Computer 
Recreations in Scientific American 
magazine) for which specimen out- 
put is provided in the form of Fiqs 5 
and 6. 

History 

This approach to pattern generation 
has been inspired by the analysis of 
the so-called Henon Attractor defined 
by (7) x r +1 = 1+y r -ax 2 r, y r +1=bx r , 
where a and b are numerical con- 
stants, r=0,1,2 . . . successively and 
starting values (x 0 ,y 0 ) are to be speci- 
fied (see, for example, M Henon, 
Communications in Mathematical 
Physics, volume 50, pages 69-77, 
1976). The chaotic situation resulting 
when a=1.4 and b=0.3 is displayed 
through four stages of magnification 
using up to 5x10 6 points with start- 
ing values (0,0). The black-on-white 
pictures presented have no artistic 
merit, but did prompt considerable 
mathematical investigations resulting 
in part from their apparent 'ultimate' 
simplicity. 

Such advantages seem unlikely in 
relation to the point generators (5) 
and (6) featured here, although some 
mathematicians are optimistic. Read- 
ers are encouraged to experiment 
with alternative forms of (4) — that 
is, to develop new designs rather 
than to seek underlying advances in 
the rather difficult mathematical 
theory. 
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